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In a recent publication, we used the data Irom WMAP and SDSS to constrain the primordial 

perturbations and to predict the B-mode polarization sourced by cosmic string networks. We have 

been alerted [jl to the existence of errors in the code U we used to calculate the Cosmic Microwave 

Background anisotropies from cosmic strings. Correcting the errors leads to a significant increase 

in the vector mode contribution to the CMB temperature and polarization anisotropies as well 

as an overall renormalization of the various string spectra. In these notes we explain the nature 

... of the errors and discuss their implications for previously published constraints on cosmic strings 

>G \ based on this code. The chief change in our results is that our derived limit for the cosmic string 

^^ . tension is strengthened: G/i < 1.8(2.7) x 10~^ at 68 (95)% confidence. We also note that the newly- 

>--^ ' enhanced vector mode contribution produces a greatly-increased amplitude for B-mode polarization 

in the CMB which could exceed the B-mode power produced by the lensing of primordial E-mode 

5-H ■ polarization into B-mode polarization. 

< 



^ 



cd 



I. INTRODUCTION 



Our collaboration has used the data from the Cosmic Microwave Background (CMB) anisotropy and from galley 

surveys to place cosmological limits on the properties of networks of cosmic strings in a series of papers |^ |^ |^ . 

We have recently learned ^ that the code we used to calculate the cosmic string-sourced CMB anisotropy and 

primordial power spectra contained several errors, the most important of which was a missing normalization factor in 

the evaluation of the vector mode power produced by cosmic strings. This latter error alone lead to a factor of ~ 8 

f~^ ' enhancement in vector-mode power. While revising the vector mode part of the code we also located an overall factor 

\^ , of 2 normalization error in all string spectra. Full details of all corrected errors are given in the Appendix. 

^D ■ The enhancement of the vector modes has implications for the results of our last paper Q, especially for the 

prediction for the B-mode polarization sourced by strings. The principal result of that work was a constraint on the 

O ' fractional power that strings can contribute to the CMB temperature anisotropy spectrum. Fortunately, this fraction 

I ' is mainly constrained by the shape of the string induced CMB spectrum, which did not change considerably as a 

r-v , result of fixing the code. Thus we believe that the bounds derived in our last paper on the parameter we called /, 

■*— * ' the fractional power sourced by cosmic strings, will remain unchanged, and can still be used to derive bounds on the 

j^ , properties of the cosmic string network. Changing the normalization of our curves, after correcting our errors, scales 

JV ' downward our fiducial string tension ^o by a factor of 1.8, so that our new fiducial string tension is G/io ~ 1-1 x 10^^. 

• 1^ ' The provisional limits on G/x that we derive using this corrected fiducial string tension are then G^ < 1.8(2.7) x 10~^ 

Si^ ' at 68 (95)% confidence. The limits we placed on cosmic string substructure, parameterized by the string wiggliness, 

?H ' ar, were weak at best j^; that determination should now be disregarded completely. We plan to perform our Markov 

Chain Monte Carlo analysis again, with the corrected string spectra, using the new three-year WMAP data |a] . Until 

we finish this analysis, we will use these provisional bounds, which we expect to be at least approximately valid 

A simple summary of what has changed in light of correcting these errors is that we now find much greater 
parity between the power in the vector and scalar modes caused by cosmic strings. Because this reapportionment of 
perturbation does not significantly alter the shape of the string sourced CMB spectrum and, hence, should not affect 
the cosmological constraints on the fraction of the total CMB power sourced by cosmic strings, the greater proportion 
of vector mode power leads to a greatly enhanced prediction for B-mode polarization in the CMB. Even with a smaller 
string tension, the B-mode power we now predict is a factor of 10 - 20 times as great as reported in our previous 
work. The factor of two uncertainty reflects our ignorance of cosmic string substructure. Given this enhancement, 
the B-mode polarization signal for a cosmologically viable network of strings could be the most prominent source 
for a B-mode signal in the CMB, with an amplitude even greater than is expected from E-to-B lensing. A host 
of experiments aimed at measuring B-modes with the relevant sensitivity (> 0.01{iiK)^ for 100 < £ < 2000) are 
currently in planning or underway [TJ. The prominence of the cosmic string B-mode prediction implies that these 
experiments provide perhaps the best opportunity for either observing a network of cosmic strings or placing limits 
on its properties. 
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FIG. 1: The scalar (blue dot), vector (red short dash), tensor (black long dash) and their sum (solid black) contributions to 
the CMB TT spectrum evaluated using the old (left) and the revised (right) codes for the cases of smooth (a,- — 1) and wiggly 
{ar = 1.9) strings with Gfi = 1.4 x 10~®. The spectra are proportional to (Gfj,)^ and this value of Gfi corresponds to the 
reference cosmic string tension used in our previous work Q, divided by \/2 (to correct an overall factor of two missing in the 
power spectra computed from the previous version of our code). Comparing the old and revised plots, it is evident that the 
new smooth string spectra have a factor of about 1.69 higher amplitude than the old spectra for the same string tension. We 
combine these two corrections - a factor of y/2 in overall normalization and a factor of « \/1.69 from enhanced vector modes 
- to infer that the correct, fiducial cosmic string tension should be G/j-o ~ 1.1 x 10~ , which is a factor of 1.8 smaller than the 
fiducial tension estimated in our previous work. The meaning of this fiducial tension, explained in our last paper, is this: G^o 
is the cosmic string tension that would be necessary for a network of cosmic strings to produce the same integrated primordial 
power as is found in the WMAP best-fit adiabatic model. The parameter a^ is the ratio of the effective mass-per-length of a 
wiggly string to that of a smooth string in the radiation era. See Refs. [a, Ig] for further details on string substructure. 



II. COSMIC MICROWAVE BACKGROUND TEMPERATURE ANISOTROPY 



As discussed in the Appendix, there was an overall factor of 2 missing in all spectra. One can correct for this 
factor by dividing the value of the string tension G/xg by a factor of V2- In particular, the value Gfio = 2 x 10~^ 
used in our previous work 5] should now be replaced by GfiQ = 1.4 x 10^^. In addition to this factor of 2, the 
enhanced vector-mode power from our corrected code leads to a larger amplitude in the TT spectrum produced by 
cosmic strings for a given cosmic string tension. In Fig. ^ we show the scalar, vector and tensor contributions to the 
string-sourced CMB TT spectrum, as well as their sum, before and after correction. We now find that vector modes 
play a much more prominent role. In the case of smooth strings their contribution is larger than that of scalar modes. 
This causes the total TT power to be nearly doubled for the same value of Gfi. Wiggliness suppresses the vector 
modes and enhances the scalar modes, a trend already pointed out in Ref. Q. The combined effect of wiggliness 
is to decrease the total TT power as compared with smooth strings. This means that constraints on Gfj, are weaker 
for wiggly strings than for smooth strings. In contrast, our previous work |a, Q reported the erroneous conclusion 
that wiggliness on strings produce an enhancement of total power and, therefore, a stronger bound on Gfi. Thus, 
the conclusions reported in those papers need to be modified. In particular, the bounds on string wiggliness reported 
there should be disregarded 

It has been accepted among the experts in the field that global strings lead to CMB anisotropics dominated by 
vector modes while local strings do not. This presumption was, to a degree, based upon the results of previous studies 
1^1^, where vector modes were found to be subdominant even for smooth strings. It now appears that this perception 
was incorrect, at least in the case of smooth strings. However, we do expect local strings to be quite wiggly, since 
they accumulate small-scale structure over the course of their evolution. Hence, we should still expect vector modes 
to be relatively low for local strings. Global strings, in contrast, are typically smooth and relativistic, so there is no 
similar mechanism available to suppress their vector modes. 
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FIG. 2: The CMB BB spectra. Left: the plot from our last paper Q showing contributions from inflation for r — 0.1 
(dash-dot), lensing (dot) and cosmic strings (solid line for smooth and dash line for wiggly). Right: revised predictions for the 
smooth (solid, red) and wiggly (dashed, blue) cases based on the assumption that strings contribute 10% of the power in the 
CMB TT spectrum. Inflation and lensing estimates are plotted as in the left hand panel but with the inflationary tensor mode 
color-coded magenta and E to B lensing still black. 

III. B-MODE POLARIZATION IN THE CMB 

Cosmic string-sourced vector mode perturbations can lead to B-mode polarization in the CMB. In view of the 
ongoing quest to observe B-type polarization in the CMB 7] , it is of particular interest to revisit the B-mode predictions 
for cosmic strings made in our previous analysis ^. In that paper we found that approximately 10% of the observed 
CMB TT spectrum can be produced by strings (fraction in strings fs = 0.1). This fraction depends only on the 
shape of the string induced spectrum, so since the shapes of the old and revised spectra, as shown in Fig.^ are quite 
similar, we will assume that fs = 0.1 is still a valid upper bound. 

Let us define B as the coefficient by which one needs to multiply the string TT spectrum in Fig. ^ in order for it to 
contribute fs = 0.1. The two spectra shown in Fig. ^ would need different values of B to satisfy this bound. In our 
previous work 5] we determined that fs — 0.1 requires i? = 0.1 in the case of (uncorrected) smooth strings (as on 
the left panel of Fig.^. We can use this information to approximate, based on the ratios between the corrected and 
uncorrected CiS, what values of B are required for the revised spectra. For the corrected smooth string spectrum we 
estimate Bgm — 0.06, while for the corrected wiggly case B^j = 0.07. 

Using these corrected B factors, we are able to produce revised estimates for the BB spectra. These are shown 
in Fig. [21 along with the plot from |^ . For a fixed G//, the net effect of correcting for the errors in the code was 
to boost the B-mode polarization from strings by a factor of approximately 20. Since the bounds on Gfi have also 
been tightened as a result of our corrections, the effective increase in the expected B-mode power is roughly a factor 
of 10 in the smooth string (a^ = 1.0) case. Because of compensating corrections, however, the enhancement factor 
remains approximately 20 in the wiggly string (a^ — 1.9) case. Both predictions appear to be well above the estimated 
limitations imposed by E to B lensing. A recent detailed study of the prospects of detection of the B-mode signal 
from cosmic strings, using a modification of our code, can be found in |13| . 



IV. CONCLUSION 



We have corrected several coding and normalization errors in the code we use to calculate the CMB anisotropy 
sourced by a network of cosmic strings. The principal changes caused by these corrections are an order of magnitude 
enhancement in the vector mode of primordial power spectra and an overall normalization correction for the string 
spectra. Tracing through the implications of these errors, we amend the principle results of our previous paper [5J 



and report a more stringent cosmological bound on the cosmic string tension: Gfj, < 1.8(2.7) x 10"'^ at 68 (95)% 
confidence. Besides tightening this bound, the increase in the vector power leads to another impHcation. Cosmic 
strings can produce B-mode polarization in the CMB much more efficiently than we estimated in our earlier work 
|5| - our previous predictions were too low by an order of magnitude. The upshot of this correction is that we now 
predict that a cosmic string network that has not been excluded by current experiments can produce a B-mode signal 
greater than any other anticipated source, including the lensing of inflationary E-mode polarization into B-mode 
polarization. This allows us to raise a tentative hope that a cosmic string network, if it exists, could be observed 
through its B-mode polarization signal. Alternatively, the non-observation of a cosmic string-sourced B-mode signal 
could provide the strictest observational limitation on the properties of such a network. 
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APPENDIX A: THE ERRORS 

1. A missing factor of 2 in the normalization of all CMB spectra 

This factor comes from the fact that the Fourier transform coefficients of the string energy-momentum tensor are 
complex numbers, while the code was evaluating only their real part. Since real and imaginary parts are statistically 
equivalent, this was OK, as along as one used the real part with the appropriate v2 correction to the amplitude. This 
factor was inadvertently omitted in previous versions of our code. 

2. A missing factor of —2\f2 in front of the vector source 



The vector part of the cosmic string CMB code [? ] follows the conventions outlined in Ref. [iJl . The conventions 
in the scalar and the tensor part of the code, written as part of CMBFAST by Seljak and Zaldarriaga, are different 
and are described in Ref. |ll| . In particular, the two conventions treat differently the derivation of the numerical 
coefficients in front of the active sources on the RHS of the Einstein equations. 

As implemented in the code, the cosmic string contribution to the vector modes comes through the Tia component 
of the Fourier transform of the energy-momentum tensor of the string network. This follows after one assumes k = zfc, 
and the full argument for why there is no loss of generality in this procedure can be found in Ref. I9l| . Here we derive 
the relation between Ti3(k, 77) (which corresponds to the variable emtV in the code) and the vector anisotropic stress 

7ri^^(k,77) that appears on the RHS of Eq. (70) of [li|. 



Ji 



In Ref. 113j ""s (k, t]) is defined as the coefficient of the expansion of T!y (x, r\) in the set of eigenfunctions Q^- (k, x) 



(1)^ 



Qg'(k,x) = 
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The relevant equation in Ref. [I^l is Eq. (40), which says 



T,, =^W(k,r;)Q«(k,x). 



zkx 



This is perhaps a bit confusing. One should, in principle, write the full expansion: 
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However, the function 



The vector part of the energy- momentum tensor T^ has contributions from both tts and 
y(k, 77) on the LHS of Eq. (70) of Ref. [10| was also defined in Eq. (37) only in terms of Q^^^'. As the authors say, 
this is acceptable, since -I- and — modes are equivalent and it is sufficient to consider just one of them, while taking 
proper care of numerical prefactors. In the code, Eq. (37) of Ref. [i3 was interpreted as 



hoj{Kij)e'^^^-VQ 



i(i) 



(A4) 



where hoj(k, rj) is the Fourier transform of hQj{'x, rf). As a result, one obtains the foUowing equation for T^(k, 77): 
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Choosing k = zfc assures that T13 and T23 are pure vector modes {i. e. they do not contribute to the scalar and tensor 
modes). For the {13} component, using Eq. (|AHl with k = zfc, we have 



Tise^^ = -Tia 
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where in the last step we used the fact that (e^, + iey)i 
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Pf^f - 2V2ri3 



1. Hence, using (|A6() in HA5|I . we can write 

(A7) 



It was the factor of —2\f2 in front of T13 that was missing in the code, as was pointed out to us by Slosar |l|. In 
other words 



7rW(k,77) = -2V2ri3(k,77). 



(A8) 



Note that this would lead to a factor of 8 increase in the vector mode contribution to the TT and BB power spectra. 
This factor would not be necessary if one followed the conventions of scalar-vector-tensor (SVT) decomposition 
used in Ref . [l^l ■ That same decomposition was used in Ref . [3 . This is probably the reason for the aforementioned 
oversight, since a cross-check against results of Ref. JJ was one of our tests. It now appears that the vector modes in 
Ref. Q were also underestimated, possibly throu gh a similar error in normalization. The tensor part of CMBFAST 
uses the SVT decomposition consistent with Ref. jl3 | and does not require additional factors in front of Ty . 



3. Other Coding Errors 

Some other errors were discovered in the vector part of the code. One was essentially a typo, whose effect was to 
set the time derivative of the scale factor during tight coupling to zero. More specifically, in subroutine f derivsv, 
instead of 



it should have been 



do 80 1 = l,nv 



do 80 1 = 3,nv . 



(A9) 



(AlO) 



The effect of fixing this error was a decrease in the polarization sourced by vector modes. 

The other errors were in subroutine f outputv, in the equations that evaluates variables dve and dvb, which are 
the E- and B-mode source functions related to variables S\ and S\ defined in Eq (61) of Ref. [iO|. The old 
equations were 

dve = (— sqrt(6.0d0)/10.0d0) * (vis(j) * polv/x * *2 + (vis(j) * polvpr + dvis(j) * polv)/(ak * x)) 

dvb = (-sqrt(6.0d0)/10.0d0) *vis(j) * polv/x (All) 

The corrected equations are 

dve = (— sqrt(6.0d0)/2.0d0) * (2.0d0 * vis(j) * polv/x * *2 + (vis(j) * polvpr -|- dvis(j) * polv)/(ak * x)) 
dvb = (-sqrt(6.0d0)/2.0d0)*vis(j) * polv/x (A12) 

Namely, there was an extra prefactor of 1/5, which was inserted thanks to a confusion caused by a difference in the 
conventions used by Refs. [lO| and [ll|, and a missing factor of 2 in front of the first term in the equation for dve, 
which was a result of a careless integration by parts. 

The combined effect of fixing the errors specified in this subsection was to increase the amplitude of both the E 
and B polarization sourced by the vector modes by a factor of about 2-3 in their power spectra. 
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